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ABSTRACT 

We present systematic spectral analyses of GRBs detected by the Fermi Gamma- Ray Burst 
Monitor (GBM) during its first two years of operation. This catalog contains two types of 
spectra extracted from 487 GRBs, and by fitting four different spectral models, this results in a 
compendium of over 3800 spectra. The models were selected based on their empirical importance 
to the spectral shape of many GRBs, and the analysis performed was devised to be as thorough 
and objective as possible. We describe in detail our procedure and criteria for the analyses, 
and present the bulk results in the form of parameter distributions. This catalog should be 
considered an official product from the Fermi GBM Science Team, and the data files containing 
the complete results are available from the High-Energy Astrophysics Science Archive Research 
Center (HEASARC). 

Subject headings: gamma rays: bursts — methods: data analysis 



1. Introduction 

Compendia of physical observations of an enig- 
matic phenomenon such as Gamma-Ray Bursts 
(GRBs) should be as complete and uniform as pos- 
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sible, so that useful inferences can be drawn from 
the assembled data. Here we have followed in a 
line of successful GRB Catalogs from the Burst 
And Transient Source Experiment (BATSE) to 
present the community with representative spec- 
tral fits for a complete sample of Fermi Gamma- 
Ray Burst Monitor (GBM) bursts, observed dur- 
ing the first two years of operation (July 14, 2008 
to July 13, 2010). The methodology of producing 
this catalog strives to balance completeness, uni- 
formity, and accuracy. The selection criteria, tech- 
niques, spectral models used for fitting, etc., all 
deriv e from a previous wo rk, based upon BATSE 
data (jMallozzi et al.lll995| ) , that is being prepared 
for publication as a complete sample from that 
earlier mission. It is our hope that the two data 
sets, the current GBM sample and the complete 
BATSE sample, can serve as a whole for compari- 
son and contrast. Further, we will continue to add 
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to this series as more observations are accumulated 
by GBM. 

Many correlative studies in the field of GRB 
global properties studies rely on only two very 
distinct representative spectra ( Amati et al]|2002t 
iGhirlanda et al.ll2004tlYonetoku et al.ll2004l ). These 
are the time-integrated and peak flux spectra, rep- 
resenting (respectively) the average emission and 
the most luminous. Defining these uniformly is 
not trivial and trade-offs inevitably must be made. 
We have chosen to accumulate every spectrum 
that meets a minimum flux significance into our 
average spectrum, which leaves out quiescent pe- 
riods that will dilute the signal. For the peak flux 
spectra we have chosen a slice of the burst repre- 
senting the peak flux on a particular time scale. 
Because of the large dynamic range of burst du- 
rations (ms to thousands of s) we have chosen 
a smaller accumulation interval for those bursts 
with durations shorter than 2 s compared to those 
longer than 2 s. In this way, we can still make a 
distinction between the fluence and the peak flux 
spectra for the traditionally short-class GRBs. 

The set of functions for fitting spectra was 
chosen for uniformity with th e existing set of 
BATSE GRB Spect r al Catalogs dBand et al.lll993t 



Kaneko et al.l l2006t iGoldstein et al" in prep. ). 



as 

well as with the (forthcoming) GBM Time Re- 
solved Spectroscopy Catalog of Bright Bursts. In 
most cases where the statistics of the spectra are 
sufficient, the data can support a number of free 
parameters that is rarely greater than four. Hence, 
we have selected spectral functions that simply ex- 
hibit two, three, and four free parameters. Our 
choice of the best model for each GRB, based upon 
the fit statistics, will be discussed below; however, 
for many cases it is clear when a simple model 
fails to fit the data, or when the fit does not sig- 
nificantly improve with a more complicated model. 

We start with a brief description of the GBM 
de tectors and ca l ibrati on i n Section 2 and refer 
to iMeeean et al.l (|2009t ) and lBissaldi et al.l (|2009l ) 
for a more thorough and complete description of 
the instrument and calibration respectively. This 
is followed, in Section 3, by a description of the 
methodology used in the production of this cata- 
log, including detector selection, data types used, 
energy selection and background fitting, and the 
source selection. We then offer a description of 
the spectral models used in this catalog in Section 



4. In Section 5, we present the spectral analysis 
methods and results. A description of the Catalog 
file format can be found in the Appendix. 

2. GBM Detectors and Calibration 

GBM is one of two instruments on-board the 
Fermi Gamma-Ray Space Telescope, which was 
launched and placed into orbit on June 11, 2008. 
GBM is a 14-detector instrument designed to 
study the gamma-ray sky in the energy band of 
—8 keV-40 MeV. Twelve of the detectors are 
sodium iodide (Nal) scintillation detectors, placed 
in groups of three on the side edges of the space- 
craft and pointed at various angles from the Fermi 
Large Area Telescope (LAT) boresight. The point- 
ing angles were optimized through simulations to 
adequately survey the entire unocculted sky at any 
time during the orbit, while also biasing cover- 
age toward the LAT pointing direction. The other 
two detectors are each composed of a bismuth ger- 
manate (BGO) crystal and are placed on either 
side of the spacecraft pointing at right angles to 
the LAT boresight. Each BGO crystal is coupled 
to two photomultiplier tubes, which collect the 
scintillation photons and convert them into elec- 
tronic signals. The Nal detectors cover an energy 
range from 8 keV-1 MeV, and the BGO detectors 
cover a range of 200 keV-40 MeV. Because of the 
all-sky coverage of the Nal detectors, and the fact 
that the diameter of each Nal crystal (12.7 cm) 
is ten times its thickness, the approximate loca- 
tion of a burst can be determined by comparing 
the relative count rates in the detectors that ob- 
served the bursts. The source location is then 
calculated in spacecraft coordinates and used in 
the production of the detector response matrices 
(DRMs). The DRMs are mathematical models of 
the detectors' response used to map the observed 
counts into photons of known energy. Each de- 
tector's response is dependent on incident photon 
energy, the measured detector output energy, and 
the detector-source angle and the earth-source- 
spacecraft geometry. 

A key issue for accurate spectroscopy is the 
fidelity of the detector calibration model, which 
matches detector pulse height analyzer (PHA) 
channels to physical energies. The individual 
GBM detectors were extensively calibrated be- 
fore launch using several nuclear line sources 
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(|Bissaldi et al.ll2009j) . [n addition to line sources, 
the BESSY accelerator in Berlin, Germany was 
used to determine the low-energy light output of 
four of the flight-qualified Nal detectors. Finally, 
a source survey calibration was performed on the 
entire assembled Fermi Observatory, just prior to 
observatory-level thermal vacuum testing. This 
last calibration consisted of several radioactive 
sources being individually placed at various an- 
gles relative to the observatory in order to verify 
and validate the spacecraft and ins trument scat- 
terin g models, used in the DRMs ( Hoover et al.l 
20081 ). 

Given that the lowest energy of the radioactive 
sources was 14.4 keV and that most of the detec- 
tors were not calibrated with a continuum energy 
source such as BESSY, it is useful to validate the 
GBM calibration using on-orbit observations of 
known astrophysical sources. A good spectroscopy 
standard is the Crab Nebula, which has been stud- 
ied since the infancy of gamma-ray astronomy. 
However, GBM is not a pointed instrument, so we 
have made use of the Earth occultation technique , 
as pioneered with BATSE ([Harmon et al.l 120021 ). 
We have obtained excellent fits that are in good 
agreement in every detector with th e canonical 
Crab Nebula spectrum (for details, see lCase et al.l 
201 lL Wilson- Hodge et al. in prep.). Notably, 
there are no systematic deviations in the residuals 
to the continuum fits that would indicate issues 
with the calibration. The Earth occultation tech- 
nique uses 16 re-binne d broad energy chan nels in 
the CSPEC data type (|Meegan et al.ll2009t ). 

When some very bright sources have been fitted 
using the 128 energy channel CSPEC data type, a 
kink in the spectrum at ^33 keV can be observed. 
This takes a 'p-cygni' type form, where some chan- 
nels (at lower energies) lie below the fitted model 
and others (at higher energies) lie above it. The 
total extent in channel space is roughly 4-6 chan- 
nels, while in energy space it is about 5-10 keV. 
Since the excess residuals have equal weight below 
and above the model fit and do not change the 
fitted parameters appreciably when the affected 
channels are excluded, we have chosen to ignore 
this effect in the current work. It is likely that the 
non-linear jump in the detector light output at 
the Na K-edge has not been model ed completely 



tion work. As an example of non-linearity of the 
response, see Figure [TJ which is the spectrum of 
the bright GRB 081009 using two Nal detectors. 

Solar flares provide an additional means of cal- 
ibration of the GBM detectors, particularly the 
high-energy BGO detectors . For example, the 
positron annihilation line (511 keV), neutron cap- 
ture line (2.2 MeV), and various nuclear lines 
(e.g., 4.4 and 6.1 MeV) were observed with BGO 
detector in the solar flare of 2010 June 12 (Abdo 
et al. in prep.), demonstrating the calibration of 
that detector. In that flare there was a small gain 
shift (1%), showing limitations of the gain adjust- 
ment software and the 511 keV line had an addi- 
tional small offset (1%), possibly showing limita- 
tions of the energy calibration of the BGO detec- 
tors. 

In addition to persistent, known sources, cal- 
ibration of GBM can be studied with other 
transient events such as Soft Gamma Repeaters 
(SGRs). SGR bursts are short, with typical dura- 
tions of several tens to a few hundred milliseconds. 
Their vF v energy spectra typically peak around 
40 keV with a very steep decline in flux above the 
peak. Analysis of all the SGR bursts detected by 
GBM is currently underway (see Lin et al. 2011; 
van der Horst et al., in preparation; von Kien- 
lin et al., in preparation). Typically the models 
and fit parameters are consistent with those found 
from the observations of other instruments (for re- 



Woods fc Thompsonl(|2^06l) : lMeregh"etti 



views, see 

(|2008l )). We fit the integrated count spectrum of 



(see t he Nal K-edge discussion in iBissaldi et al 



(2009)), which we will address in future calibra- 



one of the brightest bursts of SGR 1550 - 5418 
observed by GBM. The spectrum is best fit with a 
model that ag rees within the par ameter uncertain- 
ties with the I srael et al. (2008) results obtained 
for SGR 1900+14 with Swift). This demonstrates 
the ability of the low-energy calibration to fit the 
overall spectrum well, except for the area around 
the Nal K-edge where we notice that the cali- 
bration is not adequate in removing instrumental 
effects, as discussed above. 

Concerning the calibration of GBM with re- 
spect to GRBs, a number of efforts have been 
pursued to compare the GRB spectra attained by 
GBM to spectra obtained by other instruments 
observing the same bursts in approximately the 
same energy band. One of the first GRBs detected 
by GBM (GRB 080723B) was also observed by 
two instruments onboard the INTEGRAL space- 
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craft. A study ( von Kienlin et al. 2009) was per- 
formed to compare GBM's spectral analysis re- 
sults with the ISGR I instrument (15 keV-1 MeV; 



Lebrun et al 



200 31) and the SPI i nstrument (18 



Vedrenne et~aH 120031) . This burst 



keV-8 MeV; 
was selected because it was bright, spectrally hard, 
and seen in most GBM detectors. This analysis 
was performed before the first in-flight calibration 
was released, and before an accurate modeling of 
the spacecraft solar panels, which could block or 
attenuate the observation from some detectors. A 



later study (jTiernev et al.ll2010f ) was performed on 
four other GRBs also detected by ISGRI, taking 
into account blocked detectors and a new in-flight 
calibration. These studies have confirmed that 
GRB spectra as measured by GBM are consistent 
with other well-established and calibrated instru- 
ments. 

3. Method 

During the first two years of operation, GBM 
triggered on 491 GRBs, 487 of which are presented 
in this catalog. The remaining bursts are excluded 
due to a low accumulation of counts or a lack of 
spectral/temporal coverage. In order to provide 
the most useful analysis to the community, we 
have attempted to make the method as systematic 
and uniform as possible. When we deviate from 
uniformity we indicate the circumstances clearly. 
Here we detail the detector and data selection as 
well as the process used to fit the data. Many 
of the criteria are adopted from the GBM Burst 
Catalog, (Paciesas et al., in prep.) and we have at- 
tempted to maintain this in all aspects. However, 
due to the nature of spectral analysis we demand 
stricter criteria to ensure that we have adequate 
signal in all energy channels. This also effectively 
reduces the GRB sample from that used in the 
burst catalog. In addition, this catalog includes 
4 GRBs not included in the GBM Burst Cata- 
log. GBM trigger number 081007.224 contains an 
observati on of a GRB triggered b y Swift (GRB 
081007A; [Baumgartner et al.l2009l) at -121 s after 
GBM. The other three (GBM triggers 091013.989, 
091022.752, & 091208.623) where originally classi- 
fied as unknown triggers and were not included in 
the Burst Catalog, however their spectra are con- 
sistent with the GRB spectra observed by GBM. 
Two other triggers were classified as unknown, and 
the spectral analysis of those triggers reveal that 



they are too weak to determine if they have GRB- 
like spectra. 

3.1. Detector Selection 

GBM employs 14 detectors, situated at differ- 
ent angles to promote monitoring of the full un- 
occulted sky: twelve Nal detectors with an energy 
coverage of 8 keV to 1 MeV and two BGO de- 
tectors with an energy coverage of 200 keV to 40 
MeV. For this catalog, we use the data produced 
by both Nal and BGO detectors. GBM was de- 
signed to observe the entire unocculted sky at all 
times in order to increase the number of GRBs de- 
tected. However, when GBM triggers on a GRB, 
not all detectors have an optimal pointing for re- 
constructing the observed GRB spectrum. There- 
fore, we begin with all 12 Nal detectors and select 
a subset of these detectors with optimum viewing 
angles. To ensure sufficient detector response, we 
determined that using the detectors with viewing 
angles less than 60° from the location of a burst 
maximize the effective area and provide the best 
signal for spectral analysis, since the detector an- 
gle is strongly correlated with the count rate in the 
detector. In addition, we select the BGO detector 
that is closest to the Nal detectors in the subset. 
With this subset we run a program designed to 
determine any spacecraft blockage that would in- 
terfere with the signal. In most cases this provides 
a smaller subset of detectors that should be free 
of any blockage. However, due to small inaccura- 
cies in the detector mass model or location uncer- 
tainties, the blockage code does not always return 
a subset of detectors that is free from blockage. 
This is evident when the low-energy data deviate 
strongly from the fit model, as depicted in Fig- 
ure [2j This can affect particular detectors when 
the burst is partially occulted by the spacecraft. 
When this occurs we remove these detectors from 
the selected sample. Finally, in the event that 
more than 3 Nal detectors exist in our subsample 
for a particular burst, we discard all but the 3 de- 
tectors with the smallest source angles. This helps 
to prevent a fitting bias toward lower energies for 
bursts with more than 3 sufficient detectors, since 
the high significance of data in the Nal detectors 
will influence the spectral fitting procedure more 
than the lesser amount of data in the BGO detec- 
tors. With this final detector set for each burst, 
we have an adequate and reliable dataset for per- 
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forming spectral analysis. 

3.2. Data Types 

The primary data type used in this catalog was 
the time-tagged event (TTE) data, which provide 
absolute time resolution to 2 fxs based on GPS, 
and high energy resolution of 128 channels. For 
the purpose of this catalog, we choose a standard 
time binning of 1024 ms for bursts greater than 2 
s in duration as defined by the burst tgo (see Pa- 
ciesas et al. 2011) and 64 ms for bursts of duration 
2 s and shorter. The time history of TTE typically 
starts at ~30 s before trigger and extends to —300 
s after trigger. This timespan is adequate for the 
analysis of most GRBs. For GRBs that have ev- 
ident precursors or emissions that last more than 
300 s after trigger, we use the CSPEC data, which 
extend ~4000 s before and after the burst. This 
datatype also has a continuous 128 channel energy 
resolution which is normally accumulated at 4.096 
s temporal resolution that changes to 1.024 s reso- 
lution for —600 s starting at trigger time. CSPEC 
data was used for 22 GRBs in this catalog. 

3.3. Energy Selection and Background 
Fitting 

With the optimum subset of detectors selected, 
the best time and energy selections are chosen to 
fit the data. The available reliable energy chan- 
nels in the Nal detectors lie between —8 keV and 
— 1 MeV. This selection excludes overflow channels 
and those channels where the instrument response 
is poor and the background is high. Since the is- 
sue with the sodium K-edge noticeably affects only 
the brightest GBM bursts, which comprise a small 
percentage of this catalog, we have chosen to keep 
the channels affected by the K-edge in the analysis 
of all GRBs. The fit parameters for bursts that are 
affected by non-linearity around the K-edge do not 
noticeably change; the only evident change in the 
results is the fit statistics and hence the goodness- 
of-fit. We perform a similar selection to the BGO 
detector for each burst, selecting channels between 
-300 keV and -38 MeV. With the resulting time 
series we select enough pre- and post-burst back- 
ground to sufficiently model the background and 
fit a single energy dependent polynomial (choosing 
up to 4 th order) to the background. For each de- 
tector the time selection and polynomial order are 
varied until the \ 2 statistic map over all energy 



channels is minimized, resulting in an adequate 
background fit. 

3.4. Source Selection 

Once the background count rates are deter- 
mined for each detector for a given burst, we sum 
the count rates and background models over all 
Nal detectors to produce a single lightcurve from 
which we make objective time selections. We sub- 
tract the background, convert the rates to counts, 
and calculate the signal-to-noise ratio (SNR) for 
each time history bin. Only the time bins that 
have a SNR greater or equal to 3.5 sigma are se- 
lected as signal. The time selections produced 
from this criterion are then applied to all detec- 
tors for a given burst. This criterion ensures that 
there is adequate signal to successfully perform a 
spectral fit and constrain the parameters of the 
fit. This does however eliminate some faint bursts 
from the catalog sample (i.e., those with no time 
bins with signal above 3.5 sigma). While this 
strict cut was performed to provide an objective 
catalog, it is possible that not all signal from a 
burst was selected. However, most of the signal 
below 3.5 sigma is likely indiscernible from the 
background fluctuations, so a spectral analysis in- 
cluding those bins would likely only increase the 
uncertainty in the measurements. This selection is 
what we refer to as the "fiuence" selection, since 
it is a time-integrated selection, and the fiuence 
derived is representative of the fiuence over the 
total duration of the burst. The other selection 
performed is a 1.024 ms peak photon flux selec- 
tion for long bursts and 64 ms peak count rate 
flux selection for short bursts as defined by their 
tgo duration (Paciesas et al. in prep.). This se- 
lection is made by adding the count rates in the 
Nal detectors again and selecting the single bin 
of signal with the highest background-subtracted 
count rate. This selection is a snapshot of the 
energetics at the most intense part of the burst. 
Figure [3] shows the distribution of accumulation 
time based on the signal-to-noise selection criteria. 
The accumulation time reported is similar to the 
observed emission ti me of the burst, excludin g qui- 
escent periods, (e.g. Mitrofanov et al. 1999t ). Fig- 
ure [3] also includes the comparisons of the model 
photon fiuence and photon flux compared to the 
accumulation time. Note that both comparisons 
contain two distinct regions associated with short 
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and long GRBs. While there appears to be a very 
clear correlation between the photon fiuence and 
the accumulation time, there is little correlation 
between the burst-averaged photon flux and the 
accumulation time. However, the latter depicts a 
relationship similar to the hardness-duration rela- 
tionship explored in iKouveliotou et al.l (|1993f ) . 



4. Models 

We chose four spectral models to fit the spec- 
tra of GRBs in our selection sample. These mod- 
els include a single power law (PL), Band's GRB 
function (BAND), an exponential cut-off power- 
law (COMP), and smoothly broken power law 
(SBPL). All models are formulated in units of pho- 
ton flux with energy (E) in keV and multiplied by 
a normalization constant A (ph s _1 cm~ 2 keV" 1 ). 
Below we detail each model and its features. 

4.1. Power-Law Model 

An obvious first model choice, ubiquitous in as- 
trophysical spectra, is the single power law with 
two free parameters, 



f PL (E)=A 



E 



(1) 



where A is the amplitude and A is the spectral 
index. The pivot energy (E p i V ) normalizes the 
model to the energy range under inspection and 
helps reduce cross-correlation of other parameters. 
In all cases in this catalog, E piv is held fixed at 100 
keV. While most GRBs exhibit a spectral break 
in the GBM passband, some weak GRBs are too 
weak to adequately constrain this break in the fits 
and therefore we chose to fit these with the PL 
model. 

4.2. Band's GRB function 



Band's GRB function (jBand et al.l 119931 ) has 
become the standard for fitting GRB spectra, and 
therefore we include it in our analysis: 
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The four free parameters are the amplitude, A, the 
low and high energy spectral indices, a and /? re- 
spectively, and the vF v peak energy, E pea k. This 
function is essentially a smoothly broken power 
law with a curvature defined by its spectral in- 
dices. The low-energy index spectrum is a power 
law only asymptotically. 

4.3. Comptonized Model 

This model considered for GRB spectra is an 
exponentially cutoff power-law, which is a subset 
of the Band function in the limit that (3 — > — oo: 
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The three free parameters are the amplitude A, 
the low energy spectral index a and E pea k- E p i V 
is again fixed to 100 keV, as is the case for the 
power law model. With the extended high-energy 
response of GBM compared to BATSE, we rarely 
suffer from an inability to measure the high-energy 
spectra of medium to strong intensity GRBs; how- 
ever, we still find that the Comptonized model ad- 
equately describes many GRB spectra. 

4.4. Smoothly Broken Power-Law 

The final model that we consider in this cat- 
alog is a broken power-law characterized by one 
break with flexible curvature able to fit spectra 
with both sharp and smooth transitions between 
the low and high energy power laws. This model, 
first published in iRvdd (|l999h . where the loga- 
rithmic derivative of the photon flux is a continu- 
ous hyperbolic tang ent, has been re-parametrized 
( Kaneko et al.ll2006r ) as such: 
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In the above relations, the low- and high-energy 
power law indices are Ai and A2 respectively, E b 
is the break energy in keV, and A is the break 
scale in decades of energy. The break scale is in- 
dependent and not coupled to the power law in- 
dices as it is with the Band function, and as such 
repr esents an addit i onal degree of freedom. How- 
ever, |Kaiieko_et_alJ (|2006|) found that an appropri- 
ate value for A for GRB spectra is 0.3, therefore 
we fix A at this value. Studies on the behavior of 
A pertaining to GBM GRBs may be possible after 
more bright bursts have been detected; therefore 
this is left to possible future catalogs. 

5. Data Analysis & Results 

To study the spectra resulting from the GBM 
detectors, a method must be established to as- 
sociate the energy deposited in the detectors to 
the energy of the detected photons. This associa- 
tion is dependent on effective area and the angle 
of the detector to the incoming photons. To do 
this, detector response matrices (DRMs) are used 
to convert the photon energies into detector chan- 
nel energies. The DRM is a mathematical model 
of the deposition of photon energy in the crystal — 
a photon that interacts by the photoelectric effect 
will deposit 100% of its energy, subject to resolu- 
tion broadening, while a photon that interacts by 
a single Compton scatter may deposit only a por- 
tion of its energy. The exception to this is when a 
photon carrying the energy of the iodine K-shell 
escapes the crystal. The energy calibration de- 
termines the energy boundaries of the energy de- 
position channels. The DRMs record this energy 
response of the detectors at various angles, which 



was determined th rough extensi ve simulations us- 
ing GRESS (|Hoover et alJl2010h and validated by 
the source survey. The response matrices for all 
GRBs in the catalog were made using GBMRSP 
vl.9 of the response generator and version 2 of the 
GBM DRM database, and all responses employ at- 
mospheric response modeling to correct for possi- 
ble atmospheric interference. In particular, we use 
RSP2 files, which contain multiple DRMs based 
on the amount of slew the spacecraft experiences 
during the burst. A new DRM is calculated for ev- 
ery 2° of slew, changing the effective area of each 
detector based on its angle to the source. These 
DRMs are then all stored in a single RSP2 file 
for each detector. During the fitting process, each 
DRM is weighted by the counts fluence through 
the detector during each 2° slew segment. 

The spectral analysis of all bursts was per- 
formed using RMfit, version 3.4rcl. RMfit em- 
ploys a modified, forward-folding Levenberg- 
Marquardt algorithm for spectral fitting. The 
Castor C-Statistic, which is a modified log likeli- 
hood statistic based on the Cash parametrization 
(Cash 1979) is used in the model-fitting process as 
a figure of merit to be minimized. This statistic 
is preferable over the more traditional x 2 statistic 
minimization because of the non-Gaussian count- 
ing statistics present when dividing the energy 
spectra of GBM GRBs into 128 channels. Al- 
though it is advantageous to perform the spectral 
fitting using C-Stat, this statistics provides no es- 
timation of the goodness-of-fit, since there exists 
no standard probability distribution for likelihood 
statistics. For this reason, we also calculate \ 2 
for each spectral fit that was performed through 
minimizing C-Stat. This allows an estimation of 
the goodness-of-fit of a function to the data even 
though x 2 was n °t minimized. This also allows 
for easy comparison between nested models. 

We fit the four functions described in Section 4 
to the spectrum of each burst. The BAND and 
COMP functions are parametrized with E pea k, 
the peak in the power density spectrum, while 
the SBPL is parametrized with the break energy, 
Ebreak- We choose to fit these four different func- 
tions because the measurable spectrum of GRBs 
is dependent on intensity, as is shown in Figure 
[4] Observably less intense bursts provide less data 
to support a large number of parameters. This 
may appear obvious, but it allows us to deter- 
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mine why in many situations a particular empir- 
ical function provides a poor fit, while in other 
cases it provides an accurate fit. For example, 
the energy spectra of GRBs are normally well fit 
by two smoothly joined power laws. For particu- 
larly bright GRBs, the BAND and SBPL functions 
are typically an accurate description of the spec- 
trum, while for weaker bursts the COMP function 
is most acceptable. Bursts that have signal sig- 
nificance on the order of the background fluctua- 
tions do not have a detectable distinctive break in 
their spectrum and so the power law is the most 
acceptable function. Although for weaker GRBs 
a model with more parameters is not statistically 
preferred, it is instructive to study the parame- 
ters of even the weaker bursts. In addition, the 
actual physical GRB processes can have an effect 
on the spectra and different empirical models may 
fit certain bursts better than others. The spectral 
results, including the best fit spectral parameters 
and the photon model, are stored in files following 
the FITS standard (see Appendix for file descrip- 
tion) and are hosted as a public data archive on 
HEASARC0. 

When inspecting the distribution of the param- 
eters for the fitted models, we first define a data 
cut based on the goodness-of-fit. We require the 
X 2 statistic for the fit to be within the 3<r expected 
region for the \ 2 distribution of the given degrees 
of freedom, and we define a subset of each param- 
eter distribution of this data cut as GOOD if the 
par ameter error is with in certain limits. Follow- 
ing lKaneko et al.l (|2006l ). for the low-energy power 
law indices, we consider GOOD values to have er- 
rors less than 0.4, and for high-energy power law 
indices we consider GOOD values to have errors 
less than 1.0. For all other parameters we consider 
a GOOD value to have a relative error of 0.4 or 
better. The motivation for this is to show well- 
constrained parameter values, rather than basing 
interpretations on parameters that are poorly con- 
strained. 

In addition, we define a BEST sample where we 
compare the goodness-of-fit of all spectral models 
for each burst and select the most preferred model 
based on the difference in x 2 per degree of free- 
dom. The criterion for accepting a model with a 
single additional parameter is a change in % 2 of 



at least 6 since the probability for achieving this 
difference is ^0.01. The parameter distributions 
are then populated with the spectral parameters 
from the BEST spectral fits. 

5.1. Fluence Spectra 

The time-integrated fluence spectral distribu- 
tions admit results that are averaged over the 
duration of the observed emission. It should be 
noted that the following distributions do not take 
into account any spectral evolution that may exist 
within bursts. The low-energy indices, as shown in 
Figure [5l distribute about a —1 power law typical 
of most GRBs. Up to 33% of the GOOD low- 
energy indices violate the —2/3 synchrotron "line- 
of-dcath" , while an additional 62% of the indices 
violate the —3/2 synchrotron cooling limit. The 
high-energy indices in Figure |6] peak at a slope 
slightly steeper than —2 and have a long tail to- 
ward steeper indices. Note that the large number 
of unconstrained (or very steep) high-energy in- 
dices in the distribution of all high-energy index 
values indicates that a large number of GRBs are 
better fit by the COMP model, which is equiva- 
lent to a BAND function with a high-energy index 
of — oo. The comparison of the simple power law 
index to the low- and high-energy indices makes 
evident that the simple power law index is aver- 
aged over the break energy, resulting in a index 
that is on average steeper than the low-energy in- 
dex yet shallower than the high-energy index. We 
also show in Figure |6~(d)| thc difference between the 
time-integrated low- and high-energy spectral in- 
dices, AS = (a — (i). This quantity is useful since 
the synchrotron shock model make s predictions 



of this value in a number of cases (Preece et al 



2002). The fluence spectra distribution of AS 



bears resemblance to the time-resolved results in 
Preece et al.l ( 2002 ), a s well as time-integrated re- 



sults in Kaneko et al 



(|2006l) . 

In Figure [3 we show the distributions for the 
break energy, E ore ak and the peak of the power 
density spectrum, E peak . E break is the energy 
at which the low- and high-energy power laws 
are joined, which is not necessarily representative 
of the Ep ea k. The E oreak for the SBPL has a 
strong clustering about 100 keV, while the peak 
for the BAND distribution is closer to 200 keV. 



^ The Ep ea k distributions all generally peak around 

fottp : //heasarc . gsf c . nasa . gov/W3Bro«se/f ermi/f ermigbrst . html [ 9nn ke y an j CQver j ugt oyer twQ ord(Jrs Q f mag . 
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nitud e, which is consistent with previous find- 
ings (jMallozzi et al.lll995ULlovd et al.l|2000h from 
BATSE. As discussed in lKaneko et al.l (|2006t ). al- 
though the SBPL is parametrized with Eb re ak , the 
Epeak can be derived from the functional form. We 
have calculated the E pea k for all bursts with low- 
energy index shallower than -2 and high-energy 
index steeper than -2, and we have used the co- 
variance matrix to formally propagate and calcu- 
late the errors on the derived E pea k- The peak 
and overall distribution of E pea k is similar to that 
found by the BATSE Large Area Detectors, which 
had a much smaller bandwidth and larger collect- 
ing area. This would seem to indicate that it is 
unlikely for there to be a hidden population undis- 
covered by either instrument within the 20 keV 
- 2 MeV range. Additionally, the value of E pea k 
can strongly affect the measurement of the low- 
energy index of the spectrum, as shown in Figure 
[51 A general trend appears to show that lower 
E P eak values tend to increase the uncertainty in 
the measurement of the low-energy index, mostly 
due to the fact that a spectrum with a low E pea k 
will exhibit most of its curvature near the lower 
end of the instrument bandpass. In many cases, 
if the low-energy index is found to be reasonably 
steep (< —1), the uncertainty of the index is min- 
imized even if E pea k is low. 

It is of interest to study the difference in the 
value of E peak between the BAND and COMP 
functions since they are the two main functions 
used to study GRB spectra, and COMP is a spe- 
cial case of BAND. To study the relative devia- 
tion between the two values we calculate a statis- 
tic based on the difference between the values and 
taking into account their la errors. This statistic 
can be calculated by 



AE. 



\ E veak - E f 



peak 



peak I 



(6) 



where C and B indicate the COMP and BAND 
values respectively. This statistic has a value of 
unity when the deviation between the E peak val- 
ues exactly matches the sum of the la errors. A 
value less than one indicates the E pea k values are 
within errors, and a value greater than one indi- 
cates that the E pea k values are not within errors 
of each other. Figure [5] depicts the distribution of 
the statistic and roughly 25% of the BAND and 
COMP E pea k values are found to be outside the 



combined errors. This indicates that, although 
COMP is a special case of BAND, a significant 
fraction of the E pea k values can vary by more than 
la based on which model is chosen. 

The distributions for the time-averaged photon 
flux and energy flux are shown in Figure 1101 The 
photon flux peaks around 2 photons cm~ 2 s , and 
the energy flux peaks at 2.5 x 10~ 7 ergs cm~ 2 s _1 
in the 8-1000 keV band. When integrating over 
the full GBM spectral band, 8 keV - 40 MeV, the 
energy flux broadens and for all but the simple 
power law model, approximates a top hat func- 
tion with a small high-flux tail spanning about 2 
orders of magnitude. Note that the low-flux cut- 
off is due to the sensitivity of the instrument, and 
therefore the peak of the distributions are directly 
affected by the instrument sensitivity. A distri- 
bution of flux measurements with a more sensi- 
tive instrument will likely position the peak at a 
lower flux. Similarly in Figure II H the distribu- 
tions for the photon fluence and energy fluence are 
depicted. The plots for the photon fluence appear 
to contain evidence of the duration bimodality of 
GRBs, and have discriminant peaks at about 1 
and 20 photons cm -2 . The energy fluence in the 
8-1000 keV band peaks at about 1 x 10~ 6 erg cm~ 2 , 
while the peak in the 8 keV-40 MeV band shifts 
nearly an order of magnitude to about 1 x 10~ 5 
erg cm~ 2 for all except the COMP model. The 
COMP model is largely unaffected by the change 
in energy band due to the exponential cutoff. The 
brightest GRB contained in this catalog based on 
time-averaged photon flux is GRB 081009 with a 
flux of > 30 ph s _1 cm -2 and the burst with the 
largest average energy flux is GRB 090227B with 
an energy flux of ~ 1.3 x 10~ 5 ergs s _1 cm -2 . The 
most fluent burst (although not the longest dura- 
tion) in the catalog is GRB 090618 with a photon 
fluence of > 2600 ph enr 2 and an energy fluence 
ergs cm~^ 



of > 2.5 x 10" 4 — "™- 2 



5.2. Peak Flux Spectra 

The following peak flux spectral distributions 
have been produced by fitting the GRB spec- 
tra over the 1024 ms and 64 ms peak flux dura- 
tion of long and short bursts respectively. Note 
that the results from both long and short bursts 
are included in the following figures. The low- 
energy indices from the peak flux selections, as 
shown in Figure [TJJ also distribute about the 
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typical —1 power law. 42% of the GOOD low- 
energy indices violate the —2/3 synchrotron "line- 
of-death" , while an additional 40% of the indices 
violate the —3/2 synchrotron cooling limit, both 
of which are significantly larger percentages than 
those from the fluence spectra. The high-energy 
indices in Figure [13] peak at a slope sightly steeper 
than —2 and again have a long tail toward steeper 
indices. The number of unconstrained high-energy 
indices increases when compared to the fluence 
spectra, likely due to the poorer statistics result- 
ing from shorter integration times. As shown 
with the fluence spectra, the PL index serves as 
an average between low- and high-energy indices 
for the BAND and SBPL functions. Shown in 



Figure 13(d) is the AS* distribution for the peak 
flux spectra. This distribution is roughly con- 
sistent with the BATSE results found previously 
(|Preece et al.l 120021 iKaneko et alJl2006h . but suf- 



fers from a deficit in values due in large part to 
the inability of the data to sufficiently constrain 
the high-energy power law index. 

In Figure 1141 we show the distributions for 
Ebreak and E pea k- As was evident from the flu- 
ence spectra, the E orea k from the SBPL fits ap- 
pears to peak at 100 keV, meanwhile the E orea k 
from BAND is harder and peaks at about 200 keV. 
The E pea k distributions for all models peak around 
150 keV and cover just over two orders of mag- 
nitude, which is consistent with previo us findings 
(|Preece et alJll998tlKaneko et al.ll2006l) . It should 
be noted that the data over the short timescales in 
the peak flux spectra do not often favor the SBPL 
model, resulting in large parameter errors. Using 
the propagation of errors scheme to calculate the 
Epeak for SBPL therefore causes many of the val- 
ues to be largely unconstrained. This results in the 
much smaller E pea k distribution for SBPL seen in 
Figure |14(b)| In addition, we calculate and show 
the AE pea k statistic in Figure [15] and, as was the 
case with the fluence spectra, ^25% of the BAND 
and COMP E pea k values are found to be outside 
the combined errors. 

The distributions for the peak photon flux and 
energy flux are shown in Figure 1161 The photon 
flux peaks around 4 photons cm s , and the 
energy flux peaks at 7.4 x 10 -7 ergs cm -2 s _1 
in the 8-1000 keV band. When integrating over 
the full GBM spectral band, 8 kev-40 MeV, the 
dispersion in the energy flux increases for all but 



the simple power law model and approximates a 
top hat function with a small high-flux tail span- 
ning about 2 orders of magnitude. The GRB with 
the brightest peak photon flux is GRB 081009 at 
> 130 ph s" 1 cm" 2 . This GRB is also the bright- 
est GRB in terms of time-averaged photon flux 
which makes it an excellent burst for calibration 
studies. Similar to the time-averaged energy flux, 
the burst with highest peak energy flux is GRB 
090227B at > 8.0 x 10~ 5 ergs s" 1 cm" 2 . 

When studying the two types of spectra in this 
catalog, it is instructive to study the similarities 
and differences between the resulting parameters. 
Plotted in Figure [17] are the low-energy indices, 
high-energy indices, and E pea k energies of the peak 
flux spectra as a function of the corresponding pa- 
rameters from the fluence spectra. Most of the 
peak flux spectral parameters correlate with the 
fluence spectral parameters on the order of unity. 
There are particular regions in each plot where 
outliers exist, and these areas are indications that 
either the GRB spectrum is poorly sampled or 
there exists significant spectral evolution in the 
fluence measurement of the spectrum that skews 
the fluence spectral values. Examples of the for- 
mer case are when the low-energy index is atypi- 
cally shallow (> —0.5) or the high-energy index is 
steeper than average (< —3). An example where 
spectral evolution may skew the correlation be- 
tween the two types of spectra is apparent in the 
comparison of E pea k- Here, it is likely that a flu- 
ence spectrum covering significant spectral evolu- 
tion will produce a lower energy E pea k than is mea- 
sured when inspecting the peak flux spectrum. 

5.3. The BEST Sample 

The BEST parameter sample produces the best 
estimate of the observed properties of GRBs. By 
using model comparison, the preferred model is 
selected, and the parameters are reviewed for 
that model. The models contained herein and in 
most GRB spectral analyses are empirical mod- 
els, based only on the data received; therefore the 
data from different GRBs tend to support different 
models. Perhaps it will be possible to determine 
the physics of the emission process by investigat- 
ing the tendencies of the data to support a par- 
ticular model over others. It is this motivation, 
as well as the motivation to provide a sample that 
contains the best picture of the global properties 
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of the data, that prompts the investigation of the 
BEST sample. 

In Table Q] we present the composition of mod- 
els for the BEST samples. From this table, it is 
apparent that the fluence spectral data from GBM 
strongly favors the COMP model over the others 
in nearly half of all GRBs. The BAND and SBPL 
are favored by relatively few GRBs in the cata- 
log. It should be noted that the number of GRBs 
best fit by PL increases in the peak flux spectra 
mainly due to the fact that the smaller statistics 
from the short integration time are unable to sup- 
port a model more complex than the PL. In Fig- 
ures [18] and [19] the same error cuts used in the 
GOOD samples were also used for the BEST pa- 
rameters. Note that the PL index is statistically 
an averaging of the low- and high-energy power 
laws, and that, due to the fact that GBM spec- 
tral responses have a peak effective area at lower 
energies (below 100 keV), we have included PL in- 
dices in the BEST low-energy index distribution. 
Although in a number of cases the PL model is 
statistically preferred over the other models in this 
catalog, the spectral shape represented by the PL 
is inherently different from the shape of the other 
models. Therefore, the PL index is not necessar- 
ily representative of either the low- or high-energy 
indices from the other models. In Figures 18(a 



and 19(a) we show where the distribution of PL 



indices exists relative to the alpha and beta dis- 
tributions. The PL index represents 25% of the 
BEST fluence alpha distribution and 50% of the 
BEST peak flux alpha distribution. The fluence 
spectra, on the whole, have a steeper measured 
alpha and shallower beta than the peak flux spec- 
tra. The alpha distribution for the fluence spectra 
peaks at about —1, while the peak flux low-energy 
spectral index peaks at about —0.6. Conversely, 
the beta distribution for the fluence spectra peaks 
at —2.1 and the peak flux high-energy spectral 
index peaks at —2.4. As shown before for the 
GOOD spectra, Figure [201 shows the AS distri- 
butions for the BEST fluence and peak flux spec- 
tra. Both distributions are found to be consistent 
with previous investigations into BA TSE time- in- 



tegrated and time-resolv ed spectra (jPreece et al 



20021 iKaneko et al.|[2006h 



near 200 keV and the Eb rea k peaks at about 150 
keV, although the kurtosis of the distribution is 
much lower than that of the E pea k distribution. 
Meanwhile, the E pea k and E orea k for the peak flux 
spectra both peak at about 150 keV. As is shown 



in Figures 18(c 



Additionally, the E pea k (Figure p(b)| ) and 
Ebreak ( |19(b)| distributions differ between the 
two spectra. The fluence spectra E pea k peaks 



and 19(c) the photon flux docs 
not significantly change when widening the energy 
band from 8 keV-1 MeV to 8 keV-40 MeV. The 
fluence spectra photon flux peaks around 1.5 pho- 
tons cm~ 2 s" 1 and the peak flux photon flux peaks 
around 4 photons cm" 2 s _1 . Alternatively, the en- 
ergy flux is greatly affected by expanding the en- 
ergy band, as is shown in Figures p(d)1 and [l9(d)1 
Specifically, the distribution in the 8 keV-1 MeV 
band peak around a fewxl0~ 7 ergs cm~ 2 s -1 , but 
in the 8 keV-40 MeV band the distribution broad- 
ens significantly and becomes much less peaked, 
approaching what appears to be a top hat distri- 
bution. 

To aid in the study of the systematics of the 
parameter estimation, as well as the garner the ef- 
fect statistics has on the fitting process, we inves- 
tigate the behavior parameter values as a function 
of the photon fluence and peak photon flux for the 
fluence and peak flux BEST spectra respectively. 
These distributions are shown in Figures [2H and 
[221 When fitting the time-integrated spectrum of 
a burst, we find the low- and high-energy indices 
trend toward steeper values for exceedingly more 
fluent spectra. The simple PL index trends from 
shallow value of ~ —1.3 to a steeper value of ~ —2. 
The low-energy index for a spectrum with curva- 
ture tends to exhibit an unusually shallow value 
of ~ —0.4 for extremely low fluence spectra, and 
steepens to ~ —1.5. Similarly, the high-energy in- 
dex trends from ~ —1.6 at low fluence to ~ —2.7 at 
high fluence, although this is complicated by un- 
usually steep and poorly constrained indices that 
indicate that an exponential cutoff may result in 
a more reliable spectral fit. When inspecting the 
E pe ak as a function of photon fluence, a trend is 
much less apparent. If a burst is assumed to have 
significant spectral evolution, then obviously the 
Epeak will change values through the time history 
of the burst, typically following the traditional 
hard-to-soft energy evolution. For this reason, 
spectra that integrate over increasingly more time 
will tend to suppress the highest energy of E pea k 
within the burst, so a general decrease in E pea k is 
expected with longer integration times. However, 
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the photon fluence convolves the integration time 
with the photon flux so that an intense burst with 
a short duration may have on the order the same 
fluence as a much longer but less intense burst 
but results in a higher E pea k . This causes signifi- 
cant broadening to the decreasing trend as shown 
in Figure 21(c) The distribution of parameters 



as a function of the peak photon flux, however, is 
much less clear. The distributions shown in Figure 
|2"21 are more susceptible to uncertainty because of 
the generally smaller statistics involved in study 
the peak flux section of the GRB, except in some 
cases where the peak photon flux is on the order of 
the photon fluence. Ignoring the regions where the 
parameters are poorly constrained, another trend 
emerges from the low-energy indices; they appear 
to become slightly more shallow as the photon flux 
increases. The high-energy indices, however, ap- 
pear to be unaffected by the photon flux, except 
those that are unusually steep and indicate that 
an exponential cutoff may be preferred. Finally, 
no obvious trend is visible for E pea k as a function 
of the peak photon flux. This likely results from 
the fact that neither spectral index in the peak 
flux spectrum is affected much by the photon flux. 

5.4. x 2 Distributions 

Although the least-squares fitting process did 
not minimize x 2 as a figure of merit, we can cal- 
culate the x 2 goodness-of-fit statistic comparing 
the model to the data. To do this, we difference 
the background-subtracted count rates from the 
model rates, summing first over all energy chan- 
nels in each detector, and then over all detectors. 
This is shown by 
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where the Oij are the observed count rates, Bij 
are the background rates, My are the model rates, 
and cr 2 ^. . are the derived model variances. In the 
ideal situation and assuming acceptable spectral 
fits (i.e. when performing spectral analysis of sim- 
ulated data), the reduced x 2 value (x 2 /d.o.f.) will 
tend to distribute around a value of 1. In this way, 
X 2 gives an estimate on the acceptability of the fit. 
Figure [23] shows the reduced x 2 distributions for 
both the fluence and peak flux spectra, as well as 



the corresponding BEST distributions. It is im- 
portant to note that the distributions all peak at 
slightly larger values than 1, which is acceptable 
since the fits did not minimize x 2 and even the 
BEST spectral fits represent an approximation to 
the actual spectra of GRBs due to their empiri- 
cal nature. The fluence x 2 distributions appear 
to be much broader and shifted farther from the 
nominal reduced x 2 value than the peak flux dis- 
tributions. This is likely due to the longer time in- 
tegration intervals in the fluence spectra and the 
fact that many GRBs experience spectral evolu- 
tion. The peak flux sample captures the spectra 
of all bursts in a small slice of time at the same 
stage of the lightcurve, meanwhile the fluence sam- 
ple integrates over the duration of the emission, in 
many cases over several pulses. The proximity of 
most of the reduced x 2 values to the nominal value 
is indicative of acceptable spectral fits. 

In addition, Figure [M] plots the BEST reduced 
X 2 as a function of photon fluence and peak photon 
flux for the fluence and peak flux spectra respec- 
tively. The reduced x 2 f° r the fluence fits shows 
a marked upward trend as the photon fluence in- 
creases. The average reduced % 2 starts at approx- 
imately a value of unity at a low fluence of ^0.2 
photons cm~ 2 and increases to a value of ^2 at 
~1000 photons cm -2 . Additionally there are sev- 
eral outliers to the trend that exist at high fluence 
and exhibit even larger reduced x 2 values. This 
indicates that the goodness- of-fit is increasingly 
worse the more fluent the burst. This follows from 
the fact that in most cases extremely fluent bursts 
are long and may exhibit significant spectral evolu- 
tion, therefore the time-integrated spectral fit will 
average over the evolution and will produce a sig- 
nificantly worse fit. When inspecting the reduced 
X 2 as a function of the photon flux, we find that 
the trend is not as dramatic, with low flux spec- 
tra producing a marginally better goodness-of-fit 
than the high flux spectra. This indicates that the 
peak flux spectra do not suffer from the same sys- 
tematic effects as the fluence spectra. Most of the 
large reduced x 2 values from the peak flux spec- 
tra result from high flux bursts where the problem 
around the K-edge dominates the statistics of the 
goodness-of-fit. 
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6. Summary &; Discussion 

GBM has allowed the observation of a large 
number of GRB spectra over a broad energy range 
that has not been fully covered previously. The 
broad energy range of GBM translates to parame- 
ter distributions that are variants of distributions 



found by previous instruments (e.g. iKaneko et al 
20061 Goldstein et al. in prep, on the final 
BATSE spectral catalog). The distributions con- 
tained here are similar to those shown by previ- 
ous studies, yet contain differences that display 
the usefulness of studying GRBs over different en- 
ergy bands and sensitivities. Although the typi- 
cal average GRB spectrum has largely remained 
unchanged within the GBM energy band, the dis- 
tributions of spectral parameters are no longer as 
narrow as previously thought. We have shown in 
many cases that the fitted spectrum of a GRB de- 
pends in large part on its intensity as well as the 
detector sensitivity. This observation implies that 
weak GRBs may have the same inherent spectrum 
as their more intense counterparts, yet we are un- 
able to accurately determine their spectrum, rein- 
forcing the importance of comparing the spectral 
parameter distributions from different acceptable 
models. Additionally, the high-energy coverage 
provided by the BGO detectors allows the exten- 
sion of fine spectral analysis into a new regime. 
For example, the inclusion of high-energy data in 
many cases causes the high-energy power law to 
soften, and in some cases provides evidence for 
spectral cutoffs at even higher energies that is be- 
ing investigated by the Fermi/LAT (Omodei et 
al. in prep.). Because of the high-energy cover- 
age, we have uncovered several GRBs with peak 
energies greater than 1 MeV, extending the E pea k 
distribution to higher energies. Observations of 
high- E pea k GRBs have been observed previously 
(e.g. GRB 990123: iBriggs et al.l fl999h . but were 
typically not from time-integrated spectra. Ac- 
cording to our BEST samples, 18 GRBs have a 
time-integrated E pea k larger than 1 MeV, which 
is ^6% of the GRB fiuence spectra in the catalog 
that can be fitted with a curved function. Simi- 
larly, 8 GRBs have a peak flux E pea k greater than 
1 MeV, which is ~5% of the peak flux spectra that 
can be fit by a curved function (for details on three 
short hard bursts see iGuiriec et al]|2010l) . 

Another interesting result of the parameter dis- 



tributions is the AS parameter, the difference be- 
tween the low- and high-energy spectral indices. 
This can be an important quantity because cur- 
rent models for the GRB prompt emission mecha- 
nism caii_be_teokeninto two categories: magnetic 
(e.g. Lee et al] 120001) or in ternal /external shock 
(e.g. Rees fc Meszarosl 19921 ) driven. Implications 
about the shock driven cas e can be drawn directly 
from AS. Sar i et al ] (|l998h presented a model for 
the prompt emission spectrum where the emis- 
sion of synchrotron photons by shock accelerated 
electrons evolves due to radiative and adiabatic 
losses. A key prediction of this model is a se- 
ries of relationships between the photon spectral 
indices abov e and below the synchrotron peak. 



Preece et al.l (|2002l ) found that this relationship 



can be used to constrain the synchrotron cooling 
regime of the electrons by looking at the difference 
of the Band function spectral indices across the 
vF v peak. Therefore, investigating AS potentially 
provides insight into the emission mechanism of 
GRBs by comparing AS to the prediction offered 
by the synchrotron shock model. 

A comparison between different catalogs can 
be instructive, especially when noting the differ- 
ence in methodologies and analysis. For example, 
the method ology and fi t resul ts for GBM bursts 
presented in iNava et "aL ( 2011 ) are somewhat dif- 
ferent than the analysis and results given in this 
GBM catalog. Simple differences in data selection 
and methodology, such as analysis on detectors 
with spacecraft blockage and usage of detectors 
with detector-to-source angles larger than 60°, 
can significantly affect the res ults of the spectra l 
analysis. Additionally, in the INava et al.1 (|201lh 
sample, signal selections were not performed in a 
systematic way, but rather subjectively selected, 
based on a subjective estimate of where the signal 
fell below the background. These differences as 
well as many others can lead to large differences 
in results. We have compared the overall aver- 
age values and dis tributions b e tween this catalog 



and t he analysis of INava et al.l ([201 if ) . INava et al 



(|201ll ) report which model is most preferred for 
each burst in their sample; however they do this 
by comparing C-Stat, which does not allow for 
model comparison without Monte Carlo simula- 
tion for each burst. Unfortunately, a \ 2 goodness- 
of-fit is also not reported so there is no clue to 
the accuracy of their model selection method, nor 
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how well their models actually fit the data. This 
should also be noted as a caveat when comparing 
the two catalogs. Table [3 shows the sample mean 
and standard deviation of lNava et al.l ([201 if ) and a 



comparison to the results of this catalog, of which 
there are many notable differences. Although their 
average best low-energy index value for the time- 
integrated spectra is relatively similar, their aver- 
age high-energy index is found to be steeper and 
represents a narrower distribution, although this 
may be due to their stated inability to constrain 
the high-energy index. Their average E pea k value 
is slightly softer and does not capture as many 
of the extended high-energy E pea k values that are 
found in this catalog. The mean time-averaged 
energy flux is also softer by about 25% compared 
to our derived values. The differences in the peak 
flux spectra are even more d ivergent. The aver age 



low-energy index found by iNava et al.1 (|2011l ) is 



surprisingly shallow and the high-energy index is 
somewhat steeper and represents a much broader 
distribution than what is presented here. They 
found the average E pea k of the peak flux spectra 
to be somewhat harder, but their average peak 
energy flux is nearly an order of magnitude larger 
than reported in this catalog. These differences 
are representative of divergent methodologies and 
samples; therefore the reader should take care to 
understand how the values in the catalogs are de- 
rived. 

Additionally, when comparing ou r results to the 
BAT SE catalog of bright GRBs (jKaneko et al. 
2006), there are some marked differences as 
well, although man y of these arise because the 
Kaneko et all (2006) sample studied only the 350 
brightest BATSE GRBs. We investigated the sam- 
ple mean of the time-integrated spectra for com- 
parison (see Table [3]) and found that the average 
low-energy index was relatively the same and the 
high-energy index on average was slightly steeper. 
The difference of values in the high-energy index 
could perhaps be explained by the difference in 
spectral bandwidth between the two instruments. 
The time- averaged E pea k however, is on average 
harder, most likely due to the selection of the 
brightest bursts in that catalog. This also affects 
the Ebreak, which is consid erably harder a t over 
200 keV, and naturally the iKaneko et all {2006) 



increased sensitivity compared to GBM. 

In conclusion, we have presented a systematic 
analysis of GRBs detected by GBM during its first 
two years of operation. This catalog contains four 
basic photon model fits to each burst, using two 
different selection criteria to facilitate an accurate 
estimate of the spectral properties of these GRBs. 
We have described subsets of the full results in the 
form of data cuts based on parameter uncertainty, 
as well as employing model comparison techniques 
to select the most statistically preferred model for 
each GRB. The analysis of each GRB was per- 
formed as objectively as possible, in an attempt to 
minimize biased systematic errors inherent in sub- 
jective analysis. The methods we have described 
treat all bursts equally, and we have presented 
the ensemble observed spectral properties of GBM 
GRBs. Certainly there are avenues of investiga- 
tion that require more detailed work and analysis 
or perhaps a different methodology. This catalog 
should be treated as a starting point for future 
research on interesting bursts and ideas. As has 
been the case in previous GRB spectral catalogs, 
we hope this catalog will be of great assistance and 
importance to the search for the physical proper- 
ties of GRBs and other related studies. 
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catalog reports larger average photon and energy 
fluxes and much larger spreads in flux due to the 
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A. Spectral Catalog Data Files (SCATs) 

The spectral fit results are stored in a series of Flexible Image Transport System (FITS) data files, available 
online from NASA's High Energy Astrophysics Science Archive Research Center (HEASARCjl. This data 
product provides different types of spectral fits to be included the official spectral catalog of GBM bursts. 
The type is coded into the filename. There are two types of spectrum categories: 

1. Peak flux ('pflx') - a single spectrum over the time range of the burst's peak flux 

2. Fluence ('fine') - a single spectrum over the entire burst 

The time ranges for each category are included with the spectra. The spectra are fit with a number of 
models, whose complexity often depends on the signal-to-noise ratio of the spectrum. The current set is: 

1. Power law ('plaw'), 

2. Comptonized (exponentially attenuated power law; 'comp') 

3. Band ('band') 

4. Smoothly broken Power Law ('sbpl') 

The set of detectors used for the fit is not coded into the filename, and may not be the same as the set of detec- 
tors that triggered. The file naming convention is thus: glg_scat_all_bnyymmddf f f _type_modl_vxx . f it, 
with yymmdd = the date; f f f = fraction of a day, in thousandths; xx = version number of the file; type 
= 'pflx' or 'fine', the spectral categories; and modi = 'plaw', 'comp', 'band', or 'sbpl', the spectral model 
used. We have chosen to designate burst triggers by the fraction of a day in thousandths, because that time 
(86.4 s.) is less than the time for a complete trigger data readout, as currently defined in the GBM Flight 
Software (300 s.). Therefore, we will not confuse two successive triggers by name. 

FITS files are generally self-documenting, consisting of header and data pairs, where the ASCII headers 
are human readable and are followed by the data, as described in the header format statements, in either 
ASCII or binary format. The first (or, 'primary') header describes the file's content in general, and includes 
information that applies globally. The data section for the primary header is an optional image, and the 
GBM data files omit this. An example primary header is as follows: 

SIMPLE = T /Written by IDL: Mon Feb 14 13:36:27 2011 

BITPIX = 8 / 

NAXIS = / 

EXTEND = T /File contains extensions 

DATE = I 2011-02-14T19:36:28' / 

FILETYPE= 'SPECTRAL FITS' /Unique FITS file type name 

CREATOR = 'rmfit 3.4rcl' /Software/version creating file 

ORIGIN = 'GIOC ' /Name of organization 

TELESC0P= 1 GLAST ' /Name of mission 

INSTRUME= 1 GBM ' /Name of instrument 

0BSERVER= 'Meegan ' /Name of instrument PI 

MJDREFI = 51910.0 /MJD date of reference epoch, int part 

MJDREFF = 7.428703703703703E-04 /MJD date of reference epoch, frac part 

TIMESYS = 'TT ' /Time system 

TIMEUNIT= 's ' /Time unit used in TSTART, TST0P and TRIGTIME 

OBJECT = 'GRB081101167' / 



2 http: //heasarc .gsf c .nasa. gov/W3Browse/f ermi/f ermigbrst . html 
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DATE-OBS= '2008-ll-01T04:00:38' /Date of start of observation 

DATE-END= ' 2008-11-01T04 : 00 : 39 ' /Date of start of observation 

TSTART = 247204838.5854399 /Observation start time, relative to MJDREF 

TSTOP = 247204839.6094399 /Observation end time, relative to MJDREF 

TRIGTIME= 247204840.6334400 /Trigger time, relative to MJDREF 

FILENAME= ' glg_scat_all_bn081101167_f lnc_band_v02 . f it ' /Name of FITS file 

COMMENT This file consists of time-sequenced spectral fit parameters 

CHECKSUM= ' ZTAQbT4QZTAQbT3Q ' / HDU checksum updated 2011-02-14T19 : 36 : 28 

DATASUM = ' 0' / data unit checksum updated 2011-02-14T19 : 36 : 28 

END 

Here, we find that the event named 'GRB081101167' (OBJECT), triggered at 247204840.6334400 s Fermi 
Mission Elapsed Time (TRIGTIME), or roughly 4 hours, minutes and 38 seconds on November 1, 2008 UTC. 
In addition to this, the standard header information gives the name of the file itself, the creation time and 
program version, as well as where the file originated and details of the mission and instrument with which 
it is associated. Finally, EXTEND is true, so the primary header is followed by one or more extensions. One 
cannot assume that the order of these extensions will be preserved, however they have identifiers that will 
allow readers to access the correct one. The original GBM files will always follow the order as given here. 

The first data extension describes the detectors used and the selections made, both in time and energy 
(channel), for the spectral fit. The unique name (EXTNAME) for this extension and its structure is 'DETEC- 
TOR DATA': 



XTENSI0N= 


'BINTABLE' 




/Written by IDL: Mon Feb 14 13:36:27 2011 


BITPIX = 




8 


/ 


NAXIS 




2 


/Binary table 


NAXIS1 = 




244 


/Number of bytes per row 


NAXIS2 = 




3 


/Number of rows 


PC0UNT = 




6156 


/Random parameter count 


GC0UNT = 




1 


/Group count 


TFIELDS = 




12 


/Number of columns 


TF0RM1 = 


<20A 




/Character string 


TTYPE1 = 


' INSTRMNT' 




/Instrument name for this detector 


TF0RM2 = 


<20A 




/Character string 


TTYPE2 = 


'DETNAM' 


/Detector number; if one of several available 


TF0RM3 = 


<20A 




/Character string 


TTYPE3 = 


'DATATYPE' 




/Data type used for this analysis 


TF0RM4 = 


'20A 




/Character string 


TTYPE4 = 


'DETSTAT ' 




/Was this detector INCLUDED or OMITTED 


TF0RM5 = 


'60A 




/Character string 


TTYPE5 = 


'DATAFILE' 




/File name for this dataset 


TF0RM6 = 


'60A 




/Character string 


TTYPE6 = 


'RSPFILE ' 




/Response file name for this dataset 


TF0RM7 = 


'60A 




/Character string 


TTYPE7 = 


<FIT_INT ' 




/Fit intervals 


TF0RM8 = 


<1J 




/Integer*4 (long integer) 


TTYPE8 = 


'CHANNUM ' 




/Total number of energy channels for this detecto 


TF0RM9 = 


<2J 




/Integer*4 (long integer) 


TTYPE9 = 


'FITCHAN ' 




/Channels selected in fitting this detector 


TF0RM10 = 


= '1PE(129)' 




/ Real*4 (floating point) , variable length 


TUNIT10 = 


= 'keV 




/ 


TTYPE10 = 


= 1 E_EDGES ' 




/Energy edges for each selected detector 
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TF0RM11 = 
TUNIT11 = 
TTYPE11 = 
TF0RM12 = 
TUNIT12 = 
TTYPE12 = 
TF0RM13 = 
TUNIT13 = 
TTYPE13 = 
DATE 

EXTNAME = 

ORIGIN = 

TELESCOP= 

INSTRUME= 

OBSERVER= 

MJDREFI = 

MJDREFF = 

TIMESYS = 

TIMEUNIT= 

DATE-OBS= 

DATE-END= 

TSTART = 

TSTOP 

TRIGTIME= 

NUMFITS = 

CHECKSUM= 

DATASUM = 

END 



'1PE(128) ' 
'Photon cm~-2 s' 
'PHTCNTS ' 
'1PE(128) ' 
'Photon cm~-2 s' 
' PHTMODL ' 
'1PEC128) ' 
'Photon cm~-2 s' 
'PHTERRS ' 

'2011-02-14T19:36:28 
'DETECTOR DATA' 
'GIOC 
'GLAST 
'GBM 

'Meegan ' 



/ Real*4 (floating point) , variable length 
keV~-l' / 

/Array of photon counts data 
/ Real*4 (floating point) , variable length 
keV'-l' / 

/Array of photon model data 
/ Real*4 (floating point) , variable length 
keV~-l' / 

/Array of errors in photon counts data 
/ 



/Name of this binary table extension 
/Name of organization 
/Name of mission 
/Name of instrument 
/Name of instrument PI 
51910.0 /MJD date of reference epoch, int part 
7.428703703703703E-04 /MJD date of reference epoch, frac part 
'TT ' /Time system 

's ' /Time unit used in TSTART, TSTOP and TRIGTIME 

'2008-ll-01T04:00:38' /Date of start of observation 

'2008-ll-01T04:00:39' /Date of start of observation 

247204838.5854399 /Observation start time, relative to MJDREF 
247204839.6094399 /Observation end time, relative to MJDREF 
247204840.6334400 /Trigger time (s) relative to MJDREF 
1 /Number of spectral fits in the data 

' H7Pa J5NYH5NaH5NW ' / HDU checksum updated 2011-02-14T19 : 36 : 28 

'1048916985' / data unit checksum updated 2011-02-14T19 : 36 : 28 



Since NAXIS2 is 3, three detectors were used for this fit; the data must be read in the DETECTOR column 
to determine which ones. Descriptive ASCII text data identifies which of the GBM detectors were chosen 
(DETECTOR; e.g. 'BGCL00'), which data type was selected (DATATYPE; e.g. TTE'), the original GBM data 
file name (DATAFILE: e.g. 'glg_tte_b0_bn081101167_v01.fit'), and the time and energy selection intervals 
(FIT_INT: e.g. '-2.048: -1.024 s, 284.65: 37989. keV, channels 3: 123'). Detector-specific binary data consist 
of the PHA channel energy edges (E_EDGES), the deconvolved 'photon' counts (PHTCNTS) and uncertainties 
(PHTERRS), used to reconstruct a plot of the deconvolved fit in photon space, rather than counts space, where 
the observed data were actually fitted, using the counts model (PHTMODL). The number of energy edges is 
one more than the number of PHA channels, since they represent bin boundaries; or, physically, they are 
the energy loss thresholds in the detector, along with one more edge to complete the last bin. 

Finally, there is an extension (EXTNAME = 'FIT PARAMS') that describes the spectral fit in detail: 



XTENSI0N= 


'BINTABLE' 




/Written by IDL: Mon Feb 14 13:36:27 2011 


BITPIX = 




8 


/ 


NAXIS 




2 


/Binary table 


NAXIS1 = 




110 


/Number of bytes per row 


NAXIS2 = 




1 


/Number of rows 


PC0UNT = 




96 


/Random parameter count 


GC0UNT = 




1 


/Group count 


TFIELDS = 




15 


/Number of columns 


TF0RM1 = 


'2D 




/Real*8 (double precision) 


TTYPE1 = 


'TIMEBIN ' 




/Start and stop times relative to trigger 
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TF0RM2 = 
TTYPE2 = 
ASYMERO = 
TF0RM3 = 
TTYPE3 = 
ASYMER1 = 
TFORM4 = 
TTYPE4 = 
ASYMER2 = 
TF0RM5 = 
TTYPE5 = 
ASYMER3 = 
TF0RM6 = 
TTYPE6 = 
TF0RM7 = 
TTYPE7 = 
TF0RM8 = 
TTYPE8 = 
TF0RM9 = 
TTYPE9 = 
TFORMIO = 
TTYPEIO = 
TFORM11 = 
TTYPE11 = 
TF0RM12 = 
TTYPE12 = 
TF0RM13 = 
TTYPE13 = 
TFORM14 = 
TTYPE14 = 
TF0RM15 = 
TTYPE15 = 
TDIM15 = 
DATE 

ORIGIN = 
TELESCOP= 
INSTRUME= 
OBSERVER= 
MJDREFI = 
MJDREFF = 
TIMESYS = 
TIMEUNIT= 
OBJECT = 
DATE-OBS= 
DATE-END= 
TSTART = 
TSTOP 
TRIGTIME= 
NUMFITS = 
N_PARAM = 



'3E 

'PARAMO 
'3E 

'PARAM1 
'3E 

' PARAM2 

; 3E 
' PARAM3 



[ 2E 

'PHTFLUX 
[ 2E 

' PHTFLNC 
[ 2E 

' NRGFLUX 
[ 2E 

' NRGFLNC 
'2E 

'REDCHSQ 
'II 

'CHSQDOF 
'2E 

' PHTFLUXB 
[ 2E 

' DURFLNC 
[ 2E 

' NRGFLNCB 
[ 16E 

' COVARMAT 
'(4, 4) 
'2011-02-14119:36:28 



/Real*4 (floating point) 
/Band's GRB, Epeak: Amplitude 
/Corresponding PARAM contains (value, +/- error) 
/Real*4 (floating point) 
/Band's GRB, Epeak: Epeak 
/Corresponding PARAM contains (value, +/- error) 
/Real*4 (floating point) 
/Band's GRB, Epeak: alpha 
/Corresponding PARAM contains (value, +/- error) 
/Real*4 (floating point) 
/Band's GRB, Epeak: beta 
/Corresponding PARAM contains (value, +/- error) 
/Real*4 (floating point) 

/Photon Flux (ph/s-cm~2) std energy (8-1000) 
/Real*4 (floating point) 

/Photon Fluence (ph/cm~2) std energy (8-1000) 
/Real*4 (floating point) 

/Energy Flux (erg/s-cm~2) std energy (8-1000) 
/Real*4 (floating point) 

/Energy Fluence (erg/cm~2) std energy (8-1000) 
/Real*4 (floating point) 

/Reduced Chi~2 (1) and fitting statistic (2) 
/Integer*2 (short integer) 
/Degrees of Freedom 
/Real*4 (floating point) 

/Photon Flux (ph/s-cm~2) BATSE energy (50-300) 
/Real*4 (floating point) 

/Photon Fluence (ph-cm~2) for durations (user) 
/Real*4 (floating point) 

/Energy Fluence (erg-cm~2) BATSE energy (50-300) 
/Real*4 (floating point) 

/Covariance matrix for the fit (#free params ~2) 
/COVARMAT array dimensions 



/ 



'GI0C ' /Name of organization 

'GLAST ' /Name of mission 

'GBM ' /Name of instrument 

'Meegan ' /Name of instrument PI 

51910.0 /MJD date of reference epoch, int part 
7.428703703703703E-04 /MJD date of reference epoch, frac part 
'TT ' /Time system 

's ' /Time unit used in TSTART, TSTOP and TRIGTIME 

'GRB081101167' / 

'2008-ll-01T04:00:38' /Date of start of observation 
'2008-ll-01T04:00:39' /Date of start of observation 

247204838.5854399 /Observation start time, relative to MJDREF 
247204839.6094399 /Observation end time, relative to MJDREF 
247204840.6334400 /Trigger time (s) relative to MJDREF 
1 /Number of spectral fits in the data 
4 /Total number of fit parameters (PARAMn) 
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FLU_LOW = 10.0000 /Lower limit of flux/fluence integration (keV) 

FLU_HIGH= 1000.00 /Upper limit of flux/fluence integration (keV) 

STATISTC= 'Castor C-STAT' /Indicates merit function used for fitting 

EXTNAME = 'FIT P ARAMS' /Name of this binary table extension 

CHECKSUM= 'Vc2jVc2hVc2hVc2h' / HDU checksum updated 2011-02-14T19 : 36 : 28 
DATASUM = '3846818573' / data unit checksum updated 2011-02-14T19 : 36 : 28 

END 



The start and end times of the spectral accumulation are found in the 'TIMEBIN' array. This and the 
following columns are designed for several spectra results to be stored in a single file, as would be the case 
for time-resolved spectroscopy; for the purposes of the present work, there are only the results for a single fit. 
The following set of columns, labeled with 'PARAMn', where 'n' is an integer, contain both a fit parameter 
and the associated uncertainty. The comment for the column contains the name of the function fitted as 
well as the parameter name, separated by a colon. So, for instance, the data column 'PARAM1' contains 
the fitted Band GRB function -E pea k parameter and uncertainty. If the corresponding 'ASYMERn' keyword 
is then the 'PARAMn' columns contain symmetric errors, and if the keyword is 1 then the 'PARAMn' 
columns contain the asymmetric errors. All files created for this catalog contain symmetric errors only. 

The last parameter column ('PARAMn') is followed by nine standard analysis products, in turn followed by 
the covariance matrix for the fit. The analysis products are fluxes and fluence values, calculated over several 
different energy ranges, calculated by integrating the fitted photon model at high energy resolution. For the 
columns labelled, 'PHTFLUX', 'PHTFLNC', 'NRGFLUX' and 'NRGFLNC, the standard energy range for 
integration is given in the 'FLU_LOW and 'FLILHIGH' keyword values; in this case, 10 to 1000 keV. The 
columns with names beginning with 'NRG' indicate that the corresponding value has been multiplied by the 
photon energy before integrating. The uncertainties associated with these values have been determined by 
numerically evaluating the partial derivative with respect to each free parameter in the fit and multiplying 
these into the covariance matrix (found in the column 'COVARMAT', as a flattened row vector with n x n 
elements, n being the total number of parameters as indicated by the 'N_PARAM' and 'TDIM15' keywords. 
Rows and columns of fixed parameters (i.e.: those that have been frozen at a particular value during the 
fit) are set to zero in the matrix, as reconstructed by transforming the vector into an n x n matrix. Finally, 
the reduced 'goodness-of-fit' test statistic value is given in the 'REDCHSQ' column, as well as the reduced 
figure of merit. The reduced values of these statistics are computed by dividing the actual statistic by the 
number of degrees of freedom, stored in 'CHSQDOF'. These names are historical, since rmfit now allows 
both x 2 an d —2 log likelihood fitting. The name of the figure of merit used in the fitting process is provided 
by the 'STATISTC' keyword; which is the 'Castor C-STAT' function, in the example. 
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Fig. 1. — Spectrum of GRB 081009 using Nal 3 & 4 with a Band function fit. Note the excess in residuals 
at the K-edge. The current calibration is inadequate in removing this instrumental effect for many bright 
bursts. 
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Fig. 2. — Plot of the counts spectrum of GRB 090131 using Nal detectors 0, 9 and 10, and BGO detector 1. 
Close inspection of the ^8-60 keV range reveals that Nal shows a significant negative trend in residuals 
and a large deficit in count rate at low energy when compared to the other detectors, indicating a blockage 
from what is likely the LAT radiator, which is not properly modeled in the DRM. 
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Fig. 3. — 3(a) is the distribution of the accumulation time based on the 3.5c signal-to-noise selections. Note 
the similarity to the traditional tg distribution, with the minimum near 1 second. No other estimation of 
the duration was factored into the production of the accumulation time. |3(b)| and 3(c) show the comparison 
of the model photon fluence and photon flux to the accumulation time respectively. The fluxes and nuences 
shown in these figures are from the estimated BEST model fits. 
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Fig. 4. — Plot of the average maximum background-subtracted count rates (as a proxy for observed intensity) 
versus the number of degrees of freedom of the best fit model. The count rates in the Nal detectors for each 
burst was summed, the background subtracted, and the maximum count rate was computed. The best fit 
model was determined for each burst, and a geometric average was calculated for the maximum count rates 
of the bursts for each best fit model. The error bars shown are the la standard deviations of the distributions 
of maximum count rates for each best fit model. 
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Fig. 5. — Distributions of the low-energy spectral indices from ffuence spectral fits. 5(a) shows the distribu- 



tions of GOOD parameters and compares to the distribution of PL indices. 5(b) - 5(d) display the comparison 



between the distribution of GOOD parameters and all parameters with no data cuts. The last bin includes 
values greater than 1. 
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Fig. 6. — 6(a) - 6(c) are distributions of the high-energy spectral indices from fluence spectral fits. 



6(a) 



shows the distributions of GOOD parameters and compares to the distribution of PL indices. |6(Fj1 and|6(c) 
display the comparison between the distribution of GOOD parameters and all parameters with no data cuts. 
The first bins include values less than -6 and the last bins includes values greater than -1. |6(d)| shows the 
distribution of the difference between the low- and high-energy indices. The first bin contains values less 
than 0, indicating that the centroid value of alpha is steeper than the centroid value of beta. 
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Fig. 7. — Distributions of Eb rea k and E pea k from fluence spectral fits. 7(a) displays the comparison between 
the distribution of GOOD Eb rea k and i?b rea fe 



with no data cuts. 7(b) shows the distributions of GOOD 



E peak for BAND, SBPL, and COMP. [7(c)] and fUdj\ display the comparison between the distribution of 
GOOD parameters and all parameters with no data cuts. 



2G 




Fig. 8. — Comparison of the low-energy index and E pea k for three models from the fluence spectral fits. This 
comparison reveals a correlation between the E pea k energy and the uncertainty on the low-energy index: 
generally a lower energy E pea k tends to result in a less constrained low-energy index. 




Fig. 9. — Distribution of the AE pea k statistic for the COMP and BAND models from fluence spectral fits. 
A value less than 1 indicates the E pea k values are within errors, while a value larger than 1 indicates the 
E P eak values are not within errors. 
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Fig. 10. — Distributions of photon and energy flux from fiuence spectral fits. 10(a) and 10(c) display the 
flux distributions for the 8 keV-1 MeV band. |10(b)| and |10(d)| display the flux distributions for the 8 keV-40 
MeV band. Note that the plotted distributions contain the flux on two different timescales: 1024 ms and 64 
ms. 
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Fig. 11. — Distributions of photon and energy fluence. 11(a) and 11(c) display the fluence distributions from 
the 8 keV-1 MeV band. |1 1 (b) | and [l 1 (d) | display the fluence distributions for the 8 keV-40 MeV band. 
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Fig. 12. — Distributions of the low-energy spectral indices from peak flux spectral fits. |12(a) shows the 
distributions of GOOD parameters and compares to the distribution of PL indices. |12(b)|[l2(d) display the 
comparison between the distribution of GOOD parameters and all parameters with no data cuts. The last 
bin includes values greater than 1. 
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Fig. 13. — 13(a) - 13(c) are distributions of the high-energy spectral indices from peak flux spectral fits. 



13(a) shows the distributions of GOOD parameters and compares to the distribution of PL indices. |13(b)| 



and |13(c) display the comparison between the distribution of GOOD parameters and all parameters with 
no data cuts. The first bins include values less than -6 and the last bins includes values greater than -1. 
|13(d)| shows the difference between the low- and high-energy indices. The first bin contains values less than 
0, indicating that the centroid value of alpha is steeper than the centroid value of beta. 
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Fig. 14. — Distributions of Eb rea k and E pea k from peak flux spectral fits. 14(a) displays the comparison 



between the distribution of GOOD Et> rea k and Eb rea k with no data cuts. 14(b) shows the distributions of 
GOOD E peak for BAND, SBPL, and COMP. [14(c)] and [14(d)] display the comparison between the distribution 
of GOOD parameters and all parameters with no data cuts. 
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Fig. 15. — Distribution of the AE pea k statistic for the COMP and BAND models from peak flux spectral 
fits. A value less than 1 indicates the E pea k values are within errors, while a value larger than 1 indicates 
the Epeak values are not within errors. 
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Fig. 16. — Distributions of photon and energy flux from peak flux spectral fits. 16(a) and 16(c) display the 
flux distributions for the 8 keV-1 MeV band. |16(b)| and |16(d)| display the flux distributions for the 8 keV-40 
MeV band. Note that the plotted distributions contain the flux on two different timescales: 1024 ms and 64 
ms. 
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Fig. 17. — Peak flux spectral parameters as a function of the fluence spectral parameters. For all three 
parameters there is evidence for a strong correlation between the parameters found for the fluence spectra 



and those for the peak flux spectra. Note that the PL index is shown in both 17(a) and 17(b) for comparison 
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Fig. 18. — Distributions of the BEST spectral parameters for the fiuence spectra. 18(a) displays the selection 
of best low-energy and high-energy spectral indices. The shaded distribution depicts the location of the 
distribution of the PL index. |18(b)| shows the selection of the best i? pea fc and Ebreak- 18(c) and |18(d)| show 
the selection of the best photon flux and energy flux respectively. 
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Fig. 19. — Distributions of the BEST spectral parameters for the peak flux spectra. 19(a) displays the 



selection of best low-energy and high-energy spectral indices. The shaded distribution depicts the location 
of the distribution of the PL index. |19(b)| shows the selection of the best E pea k and Ei, re ak- 19(c) and |19(d")] 
show the selection of the best photon flux and energy flux respectively. 
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Fig. 20. — Distributions of the difference between the BEST low- and high-energy spectral indices. The first 
bin contains values less than 0, indicating that the centroid value of alpha is steeper than the centroid value 
of beta. 
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Fig. 21. — BEST fluence spectral parameters as a function of the model photon fluence. Note that the PL 
index is shown in both 21(a) and 21(b) for comparison. 
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Fig. 23. — Distributions of the reduced % 2 for each model. 23(a) and |23(b"J| are distributions of the reduced 
X 2 for each model for each burst. Note that these distributions include all fits, including unconstrained fits. 
23(c) and [23(d)] show the distribution of BEST reduced \ 2 values. The first and last bins in all distributions 
contain overflow values. Note that the peak of the distributions are slightly larger than 1. 
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Fig. 24. — 24(a) Reduced \ 2 as a function of the model photon fluence. |24(b)| Reduced \ 2 as a function of 
the model peak photon flux. 
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Table 1 
BEST GRB models 



PL 


SBPL BAND 


COMP 


Fluence Spectra 


112 (23%) 


68 (14%) 75 (15%) 


232 (48%) 


Peak Flux Spectra 


213 (44%) 


51 (10%) 69 (14%) 


154 (32%) 



Table 2 

Sample mean and standard deviation of the parameter distributions 



Model 


Low-Energy 
Index 


High-Energy 
Index 


E P ee.k(keV) 


E brea k(keV) 


Photon Flux 
(ph s -1 cm -2 ) 


Energy Flux 
(10~ 7 erg s _1 cm -2 ) 


Fluence Spectra 


PL 
COMP 

SBPL 
BAND 

BEST 


-1 54+ u ' ls 
-0 90+ 42 
i ns+ - 40 

-i-05lg:tS 


-2 29+ ' 47 

9 < 7 +0.36 
-0.47 
9 9J.+0.34 
Z - ZO -0.73 


99, 7C+483.80 
zzo. * O_ 12 3.80 

221.42+. 432 ; 04 
185.58+. 428 ? 4 
204 75+ 359 - 36 


129 74+290. 08 

1ZM - '^-74.39 

175 50+ 664 - 50 

J -'°- ou -100 .64 
199 7 ,+240.41 
J.z^. f J-_ s0 36 


2 93+3- 4 w 

2 82+ 3 - 49 
Z ' 5Z -1.34 

z - aJ -1.54 
9 09 + 3-77 

2 - 92 t?'.3l 


3.96+;^ 

4 q-,+10.54 

4 48+ 11 4 
1.88 

4.48+2V7 5 

4.03l';» 


Peak Flux Spectra 


PL 
COMP 

SBPL 
BAND 

BEST 


L -°^ — Q.2i 
n a-, +0.44 

J " uo _0.49 

-0 75+ - 41 
u - '°_0.40 
, 19 + 0.61 
1 - 1Z -0.50 


9 9 7 +0.48 
^■^ ' -0.53 

-2 16+ - 40 

zlD -0.50 
9 97+O.44 
' -0.50 


215 n3+3«>. 3 

217.88 +3 ^'72 
194 49+3i3'39 

J.3'i.'±»_ 1 p 3g 
zzo - 1 -125.92 


147 i 6 + 243 - 46 

,79 Qfi + 5086.2 
O(^.»u_ 230 02 

1 79 1 R + 253.56 
J.(^.J.Q_ 100 4g 


4 67+ 7 84 

s nq + 918 
°' Ua -2.68 

5.04+1? 9 , 
5.39± 2 ° 87 8 


5 85+ 7us 
7,29+19.52 

8.33± 2i * 7 

01+22.34 

-5.28 
oc;+22.61 





Table 3 

Comparison of the sample mean and standard deviation from different catalogs 



Dataset 


Low-Energy 
Index 


High-Energy 
Index 


fa?®) 




Photon Flux 

(ph s~ 1 cm -2 ) 


Energy Flux 
(10~ 7 erg s _1 cm -2 ) 


Fluence 


This Catalog BEST 
Nava et al. (20111 
Kaneko et al. (20061 


-1 05 + "- 44 
l.UO_ 45 

-1 02+° 49 
-1 07+ - 42 

lu 1 -0.36 


9 ,5 + 0.34 
Z - ZO -0.73 

-2.40+.° ; 24 
-2.43±g:|| 


205+^ 
260j; 23 6 


,9, + 240 
iz -J_80.4 

203+ 129 

zuo -80.0 


2.92+^ 

n 00 + 6.OI 


4.03+^ 
3.05l« 54 
8.861™ 


Peak Flux Spectra 


This Catalog BEST 
Nava et al. ("20111 


i -19 + 0.61 
- L ' 1 -0.50 

-0 55+ 42 
u - oo -0.37 


9 97+O.44 
z - z ' -0.50 
9 ,-+1.42 
z "" -0.61 


223+ 126 
24lif 3 7 


,79+254 
x ' z — 100 




°'°°-4 .98 
1 4 q + 53.8 
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